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A TREETOPS SIMULATION OF THE HUBBLE SPACE TELESCOPE 


HIGH GAIN ANTENNA INTERACTION 


Virtually any project dealing with the control of a Large Space 
Structure (LSS) will involve some level of verification by dig- 
ital computer simulation. While the Hubble Space Telescope 
might not normally be included in a discussion of LSS, it is 
presented at this workshop to highlight a recently developed 
simulation and analysis program named TREETOPS. This program 
was developed by Honeywell, Inc. under sponsorship of a 
Marshall Space Flight Center research and development program 
called Augmented Flexible Body Dynamics Analysis Program 
(AFBDAP). TREETOPS, the second program to be developed under 
AFBDAP, provides digital simulation, linearization and control 
system interaction of flexible, multibody spacecraft which 
admit to a point-connected tree topology. The HST application 
of TREETOPS is intended here to familiarize the LSS community 
with TREETOPS by presenting a user's perspective of its key 
features . 

Figure 1 outlines some of the outstanding features of TREETOPS. 
The program is intended as a tool for evaluating the interac- 
tion of a LSS and its associative control system. The program 
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FIGURE 



alleviates the control system analyst of the burden of 
generating the equation of motion (EOM) for flexible, multibody 
spacecraft. This is accomplished with the Singh-Likins formu- 
lation of Kane's method of generalized speeds. This 
formulation is amenable to the automatic generation of minimal 
order EOM which eliminate constrained degrees of freedom. The 
structure is modeled with articulated joints allowing zero to 
six degrees of freedom. Component modes are used to model 
flexible substructures, thereby permitting large relative 
angles and relative rates. Since the system equations are in 
general nonlinear, a linearization subprogram is provided which 
generates the linear system matrices suitable for control law 
development. An interactive preprocessor program is also pro- 
vided which generates and edits the input data base in an easy 
to use, menu driven fashion. 

The development history of AFBDAP is presented in figure 2. 
The first simulation program, simply named AFBDAP, was released 
in 1982. This was a proof of concept program applicable to 
structures with a chain topology. In 1984, TREETOPS was 
released featuring the linearization subprogram, an improved 
interactive preprocessor, an interactive postprocessor for 
plotting, and the ability to model structures with a tree 
topology. A limited, closed tree topology modeling capability 
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EQUALITY AND/OR INEQUALITY CONSTRAINTS 



whereby 


was provided through the use of "loop closure spring," 
spring force constraints were accommodated, but not kinematic 
constraints. The CONTOPS program was released in early 1986 as 
the final product of AFBDAP . CONTOPS admits structures with 
constrained topologies; both holonomic constraints, such as 
loop closures or prescribed velocities, and nonholonomic con- 
straints such as Coulomb damping, gimbal stops, etc. Since 
CONTOPS is still in the verification phase, this presentation 
will concentrate on TREETOPS. 

The basic TREETOPS program structure is shown in figure 3- 
Working from left to right, HITIP is the interactive preproces- 
sor used to create, edit and error-check the problem specific 
input data base, which is subsequently read into the basic 
simulation/linearization program, TREETOPS. HITIP creates and 
edits two other files, a machine dependent job control file and 
a parameter dimension file used in compiling the ANSI FORTRAN 
TREETOP program. The simulation/linearization program, which 
is usually run in a batch mode, generates three output files: 
a time-history output file, a linear system matrix file, and a 
restart file. The restart file can later be used to continue 
the simulation, perhaps to perform parameter variation studies 
about a nominal operating point. The output and linear matrice 
files can be examined with the postprocessor program TREEPLOT. 
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DATA CONTROL DISTURBANCE INTEGRATION 

LAW 














Alternatively, these files can be analyzed off-line to 
determine eigenvalues, FFT's, LQR's, etc. 

The simulation program can also interface with up to four user 
supplied files: (1) a model data input file, (2) a user sup- 
plied control law, (3) a disturbance algorithm, and (4) 
alternate integration programs as an option to the standard 
Runga-Kutta four pass. 

Figure *1 depicts the conceptual components of a TREETOPS model. 
The HITIP program allows the user to interactively select the 
components listed in figure 4 and assemble them into a 
multibody structure with control laws, sensors and actuators. 
Function generators can be treated as controller inputs or sys- 
tem disturbances. The control law can be any combination of 
continuous, discrete or user defined segments. The continuous 
and discrete controllers are composed of transfer functions (in 
S or Z), summing junctions and gain blocks. The controllers 
can be arbitrarily interconnected or tied to any of the ideal 
actuators listed. The structure is composed of flexible and/or 
rigid bodies, rotation and translational hinges, and loop clo- 
sure springs. The bodies are defined on an independent basis, 
with user defined local coordinate frame locations and orienta- 
tions. The hinges are used to define Euler angles for 
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FIGURE 



coordinate transformations. 


Like the actuators, the sensors 


also represent ideal devices while the sensor and actuator 
dynamics can be built into the control law definitions. 

Figure 5 presents a segment of the input data base generated by 
HITIP for the HST TREETOP simulation. Body 1 data is identi- 
fied as a flexible body which has seven node points and six 
modes. Hinge 1 is a "fictitious hinge" which connects the base 
body of a model (always labeled at 1) to the inertial frame, 
herewith six degree of freedom. Hinge 2 is a single degree of 
freedom hinge which connects body 2 to body 1 . All of the 
input data appears in a similar format, wherein the user merely 
fills in the blank of the menu driven preprocessor. Table 1 
describes the set of modal data required for each flexible 
body . 

Figure 6 depicts the TREETOPS model of the HST with two High 
Gain Antennas (HGA's) mounted on flexible masts. Body 1 repre- 
sents the core of the HST along with the two solar array 
panels. Bodies 2 and 5 represent 3.27 meter, deployable masts 
which also serve as waveguides for the antenna RF signals. 
Each antenna is connected to the mast through a pair of orthog- 
onal, single axis gimbals with dc torque motors. The table on 
the lower right identifies the vibration frequencies of each 
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Table 1 . Modal Data for the j^ 1 Body 


Symbol 

Mnemonic 

Dimension 

Description 

j 

IDBODY 


Body ID number (j^ body) 

NM 

NMODE 


Number of modes for body 

NN 

NNODE 


Number of nodes for body 

4>^M<t> 

MMASS (k, k) 

(NM, NM) 

Augmented modal mass (kg) 

4> t D4> 

MDAMP (k, k) 

(NM, NM) 

Augmented modal damping (nt/msec) 


MSTIFF (k, k) 

(NM, NM) 

Augmented modal stiffness (nt/m) 

a-'k 

ALPHABJ (i, k) 

(3, NM) 

Mass center mode shape of k““ mode 

*>k to?) 

PHIDNLBJ (i, k J) 

(3, NM, NN) 

Mode shape of k** 1 2 3 mode at the/** 1 node 

(rf> 

PHIPNLBJ (i, k ,/) 

(3, NM, NN) 

Mode slope of k^ mode at thetf^ node (rad/ 
m) 

Wk 

RJXPHIBJ (i, k) 

(3, NM) 

h-Parameter due to the k^ mode (kg-m) 

Mi,* 

ss k 

DINERTM (ICOL, IROW, k) 

(3, 3, NM) 

Change in inertia of body due to the k^ 
mode (kg-m) 

fV 

PDYADIC (ICOL, IROW, q, k) 

(3, 3, NM, NM) 

Change in inertia of the k^ mode due to 
the q® 1 mode (kg) 

-V 

PHIXPHIJ (i, q, k) 

(3, NM, NM) 

Coupling influence of the q^ mode into the 
k™ modal equation (kg) 


NOTES 


1. Index Definition * k = Mode, i = Axis, J ? = Node, q = Mode, IROW, ICOL = Dyadic matrix elements 

2. * indicates optional data is deleted when “time varying inertial option = no.” 

3. ** indicates optional data is deleted when “PHIXPHI modal term option = no.” 
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HST/HGA TREETOPS MODEL 
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HGA and the core HST. The basic problem to be addressed is the 
interaction of the fundamental antenna mast modes at 0.531 with 
the 0.515 HST mode. The stringent 0.007 arc-second pointing 
stability requirement of the HST budgets only 0.003 arc-second 
disturbances of the Line of Sight (LOS) to be caused by the 
HGA's. Analysis showed the possibility of Dahl friction in the 
HGA gimbal bearings to excite the fundament mast modes during 
slow HGA tracking maneuvers. TREET0PS was employed to verify 
the predictions made by the HST prime contractor. 

Figure 7 shows the detailed components of the HST TREET0PS 
model, consisting of seven bodies three of which are flexible, 
seven hinges with a total of three translational and nine rota- 
tional degrees of freedom; eleven sensors, seven actuators, a 
discrete and a user controller. The discrete controller models 
the HST fine guidance control system. The user controller mod- 
els the HGA discrete P-I-D controllers along with bearing 
friction and roughness, cogging torques and quantization non- 
linearities. The model thus contains 22 degrees of freedom 
plus 14 controller states. 

Figures 8, 9, and 10 show results obtained by TREET0PS. In 
figure 8, a TREET0PS HGA step response compares favorably with 
simulation results obtained by the HST contractor. This plot 
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HST/HGA TREETOPS SIMULATION 
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shows the net torque applied between bodies 2 and 3 of figure 6 
in response to a two degree step command (note the difference 
in torque units). Figure 9 shows the time history of the HST 
LOS during a typical HGA tracking maneuver of the Tracking and 
Data Relay Satellite. Notice that the 0.003 arc-second error 
budget is exceeded at the end of the maneuver as the HGA point- 
ing control system overcomes the Dahl friction forces at very 
slow speeds. The 480 seconds of simulation time required 
approximately 56 minutes of CPU time on a Sigma V computer. 
Figure 10 summarizes the LOS disturbance caused by the HGA’s 
over the same maneuver. The interaction of the HGA mast modes 
with the HST is again demonstrated as the error budget is 
exceeded. Note also that the data plotted in figure 8 is stan- 
dard TREETOPS output data while the information in figures 9 
and 10 were generated during run time in the user controller 
subroutine, which increased CPU time somewhat. 

The proceeding discussion indicates some of the potential 
applications of the TREETOPS program. Typically, programs of 
the magnitude receive a rather cool response until a sufficient 
number of users and check cases have been used to eliminate the 
unavoidable software bugs and idiosyncrasies. At MSFC, 
TREETOPS has been employed on the six programs listed in figure 
11. In each case alternate simulation results were available 
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FIGURE 11 


for verification. Honeywell also has an extensive list of 
check cases and projects used to verify TREETOPS. Over twenty 
other companies, universities, or other Government agencies 
have received a copy of TREETOPS and are currently evaluating 
its potential applications. Most responses indicate close 
agreement with other such programs while acknowledging 
TREETOPS' easy to use format and versatility. Soon CONTOPS 
will be released to the public domain through NASA's COSMIC 
network whereby maintenance provision should become available. 
Figure 11 also list some future enhancement considered for 
TREETOPS and CONTOPS. 
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